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In this paper, starting from a lattice model of topological insulators, we study 
the quantum phase transitions among different quantum states, including quantum 
spin Hall state, quantum anomalous Hall state and normal band insulator state by 
calculating their topological properties (edge states, quantized spin Hall conductiv- 
ities and the number of zero mode on a 7r-flux). We find that there exist universal 
features for the topological quantum phase transitions (TQPTs) in different cases 
: the emergence of nodal fermions at high symmetry points, the non-analytic third 
derivative of ground state energy and the jumps of the topological "order parame- 
ters". In particular, the relation between TQPTs and symmetries of the systems are 
explored : different TQPTs are protected by different (global) symmetries and then 
described by different topological "order parameters". 



I. INTRODUCTION 



As the first topologically ordered phase, the Quantum Hall Effect (QHE) is a remarkable 
achievement in condensed matter physics 1 ^. In QHE state, at low temperatures and in 
strong magnetic fields, quantized Hall conductance can be observed due to the Landau levels 
formed by states of a two-dimensional electron gas. Just one year later, people discovered the 
so-called anomalous Hall effect where additional Hall current in the ferromagnetic material 
was observed^. Then, people find that the transverse transport also exists for different spin 
species, which is called as spin Hall effect^. Recently, the topological insulators including 
the quantized spin Hall (QSH) effect^- and the quantized anomalous Hall (QAH) effect^^ 
become hot issues. 

On the one hand, to describe the quantized anomalous Hall state, a class of topologi- 
cal insulators with time-reversal symmetry breaking, the Chern number, so called TKNN 
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integer, is introduced as integrals over the Brillouin zone (BZ) of the Berry field strength™ 

Q = --!-/ (V fc x A) z dk x dk y . (1) 
^ Jbz 

Here Aj(k) = —i(uj\Vk\uj) is a Berry connection for single-electron Hamiltonian with the 
periodic part u(k) of a Bloch state ?/>k = w(fc)e* k ' r . 

On the other hand, for quantized spin Hall state the spin Chern number Q s = — Q± 
is proposed in Ref.— to be the topological invariant to characterize the topological insu- 
lators. However, due to spin mixing term (the Rashba term), the spin Chern number Q s 
is not well defined. Thus for a class of topological insulators with time-reversal symme- 
try (T-symmetry), due to Kramers degeneracy, Kane and Mele proposed a Z 2 topological 
invariant^ 1 — ~— ■ 

4 N 

(- i ) A =nn^(^), (2) 

i=l m=l 

— * — * — * 

where ki are the four high-symmetry points satisfying k = (0, 0), (0, ir), (?r, 0), (ir, tt), £,2m{ki) 
is the parity eigenvalue at each of these points, and N is the number of Kramers pairs below 
Fermi surface. Such Z 2 topological invariant can be also defined in terms of the spin Chern 
number 

A = — I d 2 k (V fc x A) 2 mod2. (3) 

2 71 " J EBZ 

Here the integral of the Berry field strength is defined on half of the Brillouin zone. Physi- 
cally, A is identical with the number of pairs of helical edge modes. 

Thus, an interesting issue is the nature of the quantum phase transitions between different 
types of topological insulators. It is well known that in Landau theory different orders 
are classified by symmetries. The phase transitions accompanied with (global) symmetry 
breaking are always described by order parameters. However, Landau's theory fails to 
describe the topological quantum phase transition (TQPT). Such type of quantum phase 
transition cannot be classified by symmetries^. Instead, they may be characterized by some 
topological "order parameters", such as the Chern number or Z 2 topological invariant^ - —. 

In this paper, we study the TQPTs between topological insulators and show their physical 
properties. We find that there are universal features of the TQPTs for different cases : the 
existence of nodal fermions at high symmetry points, the non-analytic third derivative of 
ground state energy and the jumps of the topological "order parameters". In particular, 
we find the symmetry protected nature of the TQPTs : different TQPTs are protected by 
different (global) symmetries and then described by different topological "order parameters" . 
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The remainder of the paper is organized as follows. In Sec. II, the two-dimensional lattice 
models of topological insulators are given. In Sec. Ill, we focus on the TQPTs between 
different quantum states with ^^-conservation (without time reversal symmetry). In this 
section, the global phase diagram and the critical behavior are obtained. In addition, the 
topological properties of different quantum states are calculated, including the edge states, 
the quantized Hall conductivities and the induced quantum numbers on a tt-Aux. In Sec. IV, 
we will study the TQPTs with time reversal symmetry (without S^-conservation). In Sec.V, 
We will study the TQPTs without time reversal symmetry and S^-conservation. Finally, 
the conclusions are given in Sec. VI. 



II. THE LATTICE MODELS OF TOPOLOGICAL INSULATORS 

As a starting point, we consider a lattice Hamiltonian for two-flavor spin-1/2 fermions^: 

H 1 = H + H D , (4) 



where 




fi/2 
fi/2 




•qj + X^ r z - \l/ t 

*> CT ~ / 4 'i,i+x * i+x,a 

%,a + H.C., (5) 



(6) 



—t isgn(a)/2 
2sgn(a)/2 t 

-sgn(cr) 
sgn (a) 

Here ty\ a = (V'Iat; V4b<t) * s a two-component particle operator where a =t, 4 is the spin 
index, A, B denote the flavor indices and i labels the site on the square lattice, t is the 
real hopping parameters. \i is the orbital splitting energy. D is an "effective" magnetic 
field which was first introduced in Ref.-^i phenomenologically which may be induced by local 
magnetic moments. The bond variables rf • is set to +1 everywhere at the beginning. For 
the case of D > 0, we set the sgn(cr)-D is positive for spin-up while the sgn(a)D is negative 
for spin down. Without D terms, the model in Eq.[3] is really a lattice realization of the 
Kane-Mele model proposed in Ref.— of which the ground state is a topological insulator 
for t > 1/4. 



Using the Fourier transform = -^Yl l e lk Ri ^k: w e can transform the four field opera- 

k 

tors on each site into momentum space ^ = At , \&£ Bt , ^j, \P£ B ^). Then we get 



/i_ 



where, for the spin up part 



+(*) 



— j — 2t(cos k x + cos fcj,) — D sinky + i sin 



sin k y — i sin ^ + 2t(cos /c x + cos k y ) + -D 

and for the spin down part 



h. 



(fc) 



— | — 2t(cos A^ + cos k y ) + D —(sin k y — i sin k x ) 

— (sin fcj, + i sin k x ) | + 2t(cos A^ + cos — £> 

The energy spectrum is then given by diagonalizing the Hamiltonian 

E k = (sin 2 k y + sin 2 k x ) + + 2t(cos k x + cos ± D) 2 . 

It is obvious that the term of sgn(cr)D breaks the T-symmetry, 



(8) 



H^ty^Q^HK-tyQ (9) 

where 6 = ia y K (<J Xt y !Z are Pauli matrices and K stands for complex conjugation). However, 
the model has an additional spin rotation symmetry around z-direction as 

H 1 = e~ i0<T 'H 1 e ie<T '. (10) 

In addition, by adding the Rashba term to H , we get another lattice model H 2 = H +H R 
where 

H R = R(*U*»*A ~ + E iR K^^ - + H.c. (11) 

i i 

By using the Fourier transform we get 

^-^ + / i?(sin k v + i sin k x )r° \ 

k \ R(sin k y — i sin k x )r / 

where, r° is the unit matrix. Now there is T-symmetry 

H 2 (k) = e- 1 f/f(-ife)e. (13) 
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However S z is not a good quantum number due to 

H 2 ^e- iea *H 2 e iea *. (14) 

Furthermore, one may consider a lattice model H 3 with both H D and H R as if 3 = 
ifo + Hd + iJjj. Now the Hamiltonian H3 has neither T-symmetry nor ^-conservation due 
to 

H 3 (k) ^ Q-'Hji-k)® (15) 

and 

H 3 ^e- i6a *H 2 e iea *. (16) 

In the following parts we will study TQPTs by keeping different symmetries based on 
different lattice models #1,2,3 • 



III. TQPTS WITH ^-CONSERVATION : D + 0, R = 
A. Global phase diagram 

In this section, to learn the TQPT of topological insulators with S^-conservation (D^O, 
R = 0), we focus on H\ = H + Ho. 

By calculating the phase boundary with zero fermion energy, we obtain the phase bound- 
ary, 

l±4t = ±D. (17) 

In the global phase diagram which was shown in Fig.l, there exist three phases : QAH 
state, QSH State and NI state. Taking D = 0.5 as an example (the horizontal dash dot line 
in Fig.l), one can see that the QSH state occurs in the region of t > 0.375. Although the 
T-symmetry is broken in this region, the topological properties of the QSH state is preserved 
(the edge states, the quantized Hall conductivities and the induced quantum numbers on 
a 7r-flux). With the decrease of t, the system turns into a QAH state in the region of 
0.375 > t > 0.125. In this region, one may find that the topological properties are different 
from that in the QSH state. In the region of 0.125 > t > 0, the ground state turns into a 
normal band insulator with trivial topological properties. In Fig.l, the solid, dashed and 
dotted lines are given by — 1 + At = D, 1 — At = D and 1 + At = D which denote the phase 
boundaries obtained above. 
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FIG. 1: The global quantum phase transition diagram. The solid, dashed and dotted lines are 
— /u/2 + 4t = D, n/2 — 4t = D, /x/2 + 4i = D respectively, these lines consist of the boundary of the 
quantum anomalous Hall state. The solid line is the critical line separating the QAH and QSH. 
The dashed and dotted lines are the critical lines separating the QAH and NI. The dash dot line 
intersects the solid and dotted lines at the points (t, D) = (0.125, 0.5) and (t, D) = (0.375, 0.5). 
The parameter [i is set to be 2. 

B. Topological "order parameter" - Spin Chern number 

When D ^ 0, one can not use the Z 2 topological number to classify the topological 
insulators. However S z is still a good quantum number, the Hamiltonian is decoupled 
for each spin component, then the spin Chern number is suitable to characterize different 
quantum states. The Chern number of each spin component is defined as Q^,i- 

Now we write the Hamiltonian of up-spin and down-spin fermions into 

H = E *lt N t(*o*^ + E *I n ^)*h, (is) 
k k 

where \l/fc t = (^k,At, ^k,Bt) T and ^ H = (^k,Ah ^k,Bi) T - In above Hamiltonian, we have 
define 

N t (fc) = N^t* + N^ 2 r y + A^r 2 , (19) 
N±(fc) = Nl^ + N^ry + N^, 
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where 

jV ti 3 = -1 - 2t[cos k x + cos ky] - D, (20) 
-W"t.2 = — sinfe^, N^x—sinky, 

and 

A^ >3 = -1 - 2t[cos k x + cos fey] + D, (21) 
iVj.,2 = -sinfe^, A^,i = - sin k y . 

r x ' y,z are Pauli matrices. 

Now two Chern numbers and Q± for up-spin particle and down-spin particle are 
obtained as^, 

<2t = A / rf2A; [ n t • ^ n t x <V t ], (22) 

and 

^ = A / • 9 * n 4- x ^ n J- (23) 

71 " Jn 

Here n-f- and are defined as n-f- = and = p^, respectively, is the area of Brillouin 
zone. 

In the QSH state, for up-spin, the Chern number is = 1; for down-spin, it is Q± — — 1. 
So the total Chern number is zero Q = + Q± = 0. However there is nonzero spin Chern 
number Q s = — Q± = 2. In the QAH state, for up-spin, the Chern number is — 1; 
for down-spin, it is Q4. = 0. So the total Chern number is obtained as Q = + Q± = 1 
and spin Chern number Q s — — — 1. In NI state, the Chern number is zero for 
each spin component, Q = = Q± = 0. From the results of the Chern number of different 
quantum states, one may imagine that the QAH state may be considered as half of the QSH 
state. This statement is also supported by the following calculations of other topological 
parameters. 

In order to give a clear comparison of the total Chern number and the spin Chern number 
in different quantum states, we list the results in Table. 1. 





QSH 


QAH 


NI 


Total Chern number 





1 





Spin Chern number 


2 


1 
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From the Table. 1, one can see obviously that the spin Chern number can be regarded as 
the topological invariant to characterize the topological insulators and their quantum phase 
transitions for the Hamiltonian with S^-conservation 13,25 . 



1. Edge states 

In this part, we calculate the edge states in different states. 

From the relationship between the spin Chern number and the number of the (chiral) 
edge states for each spin component n^, we have 

n t = — — [ d 2 k n t • d x ii t x d y n t , (24) 
m Jn 

n i = d2k n i ■ d x n l x d y n i- 

Different Chern numbers also characterize different edge states. In the QSH state, the total 
number of the edge states is obtained as^- 

n = n-f- + rii = 0, (25) 
n s = n t - = 2. 

Thus there are two edge states with opposite chirality. In the QAH state, there is only one 
edge state n — + 714, = 1. The number of edge states in different phases are plotted in 
Table.2. 





QSH 


QAH 


NI 


Edge state's number 


2 


I 






The numerical results of the edge states in different states are shown in Fig.2. From the 
Table.2, one can see that it is the spin Chern number that determines the number of the 
edge states in different phases. 



2. Quantized Hall conductivity 
Next, we study the quantized Hall conductivity of different states. 

We could calculate the charge Hall conductance a xy and the spin-Hall conductance a StXy 
for this model in different phases using the standard representation 21 

&xy = aJ t,xy + a i,xy (26) 
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FIG. 2: (A): Energy spectrum of the QSH state when the open boundary condition is imposed in 
x-direction. The parameters are D = 0, t = 0.75, \i = 2. (B): Energy spectrum of the QAH state 
when the open boundary condition is imposed in x-direction. The parameters are D = 2, t = 0.5, 
H = 2. 



and 



where 



&s,xy &~[,xy &\.,xy. 



d 2 k ■ d x ii^ x d y ri^, 



(27) 



(28) 



n 



/ d 2 k n 4 • d x n± x c\n 4 . 
Jn 



' xy 8tt 2 /i 

Here one can see that the physical consequence of the edge states is just the quantum 
Hall conductance in our case. In the QSH state, the quantum Hall conductance is obtained 
as 

(29) 



a t,xy — ^ , °~i,xy 



e 

T 



Thus there is no charge Hall conductance, a xy = u^ <X y + a i,xy = 0. Instead, the spin-Hall 
conductance is non-zero, cr s ,xy = v^xy — a i,xy = ^f-- In the QAH state, the quantum Hall 
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conductance is obtained as 

e 2 

Thus both charge Hall conductance and spin-Hall conductance in QAH state are half of the 
spin-Hall conductance in QSH state 

@xy — &"l,xy ' ^i,xy — i (~-U 
&s,xy — a t,xy a l,xy — ^ • 

It is pointed out that the spin-Hall conductance is still well defined due to ^-conservation 
in the QAH state. 

The spin Hall conductance in different states were shown in Fig. 3. The x-axis corre- 
sponded to the horizontal dash dot line in Fig.l. We can see that the spin Hall conductance 
is 0, and ^ in NI, QAH and QSH respectively. The spin Hall conductance changed from 
to v at the phase transition point (t = 0.125) and changed from ^ to ^ at the phase 
transition point (t = 0.375). In this sense, the TQPTs between QSH, QAH and NI are 
similar to the quantum Hall plateau transition from Dirac fermion in Ref.— . 



CD 



D 1 



0.00 0.125 0.25 0.375 0.50 



0.75 



1.00 



FIG. 3: The spin Hall conductance with respected to the three phases: NI, QAH and QSH. The 
x-axis corresponded to the horizontal dash dot line in Fig.l. 

The clearly comparison of charge-Hall conductance and spin-Hall conductance in different 
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states are shown in Table. 3 





QSH 


QAH 


NI 


Charge Hall conductance (4r) 





1 





Spin Hall conductance (4-) 


2 


1 






This table indicates that spin Chern number and total Chern number determine charge Hall 
conductance and the spin-Hall conductance in different states, respectively. 

3. Induced quantum numbers on a ir-flux 

In addition, we calculate the induced quantum number on a tt-Aux in different quantum 
states. A 7r-ffux denotes half a flux quantum (|<&o> = ) on one plaquette of the square 
lattice which is shown in Fig.4. 





E 
































-1 








i — i 


-1 


i i 





-1 



FIG. 4: Illustration of a 7r-flux in a two-dimensional square lattice ( the gray area ) at the end of 
a string of minus signs ( — 1) of links. 

We calculate the induced quantum number on a 7r-flux in the QSH state firstly. In the 
QSH state, the 7r-flux defects defined above may be considered as the Jackiw-Rebbi soliton^Z. 
From the numerical results, we find that there are two zero modes of a single 7r-flux defect. So 
such a two dimensional point defect possesses four soliton states, two for each spin (occupied 
and non-occupied) which are denoted by 

|t+>®|++>> It-)® II-), (32) 

and 

\t-)®\U) |t+)®l+->- (33) 
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Here and ||_) are the empty states of the zero modes \l/°(r) and \??(r). Around a it 
flux, the fermionic operators are expanded as 

if a (r,t) = J2^e- lEkt Vak(r) (34) 

where 6 a fc and d\ k are operators of k ^ modes that are irrelevant to the soliton states 
discussed below. a° a are annihilation operators of zero modes and a is the spin index. Thus 
we have the relationships as 

a? | t+)=lt->> a? It-) = 0, (35) 

a; I i+)=\i-), a; II-) = o. 

We define the induced fermion number operators of the soliton states, N^p with 
N a ,F = (a°Ja° a + ^Uah - $ ak d ak ) - \ 

From the relation in Eq.f l35|) . we find that |t±) or |^±) have eigenvalues of ±| of the induced 
fermion number operators, 

%A t ±) = ±~| t±>, #mH l±> = 0, (36) 

Ni,f\ I ±) = ±^\i±), N i>F \ t±) = 0. 

From N at F, we can define two induced quantum number operators, the total induced fermion 
number operator Np = N^f + N^f and the induced quantum spin operator, S z = ^(N^ t p — 

Then we calculate two induced quantum numbers defined above. For an insulator state, 
the ground states are the two degenerate soliton states denoted by |t-)® |4-+) and |t+)<8> ||- 
). One can easily check that the total induced fermion number on the solitons is zero from 
the cancelation effect between different spin components 

Nf It-)® 11+) = N F |t+)® II-) = 0. (37) 
On the other hand, there exists a spin-j moment on the soliton states |t-)® ||+) and 
lt+>® II-), 

§' I t-)® 11+) = \ It-)® II+), (38) 
S z | t+>® II-) = ~\ lt+>® II-)- 
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The induced spin moment may be straightforwardly obtained by combining the definition 
of S z and Eq. fl36|) together. 

The four different ways of occupation these two zero modes give rise to four different type 
of fluxons with the following quantum numbers: (charge = 0, spin — ±|) and (charge = ±1, 
spin = 0). The results is consistent to that there exist spin-charge separated solitons in the 
presence of 7r-flux with induced quantum numbers^ 1 ^ 1 ^. 

The charge density of a fluxon in the QSH state on a 24 x 24 lattice ^(r)^(r) + 
\I/°(r ) (r ) is shown in Fig. 5. The zero mode is localized around the defect center within a 
length-scale ~ m -1 . Here m is the mass gap of the fermions. 




FIG. 5: The charge density of a fluxon in the QSH state on a 24 x 24 lattice with periodic boundary 
condition. The charge bound to the defect is ±e. 

Secondly, we study the zero modes and the induced quantum numbers on a 7r-flux in the 
QAH state. In the QAH state, there is only one zero mode of a single 7r-flux defect. The 
charge densities \l/°(r)\l/°(r) + \I/°(r)\[ / °(r) of a pair of defects in the QAH state on a 24 x 24 
lattice with periodic boundary condition was shown in Fig.6. For the two defects, the two 
zero modes slightly split due to tunneling effect between them. 

Thus there are two soliton states, two for up spin particles (occupied and non-occupied) 
which are denoted by \\ + ) and |t-)- I n contrast, there is no zero mode and soliton states 
of the 7r-flux defect for down spin particles. Similarly we can obtain the eigenvalues of the 
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FIG. 6: The charge densities of a pair of defects in the QAH state on a 24 x 24 lattice with periodic 
boundary condition. The charge bound to the defect is ±| . 

total induced fermion number operator N F , 

Nf\U) = N^f\U) = ±\\W- (39) 
And the induced quantum spin number on the soliton states \f±) is obtained as 

S z I t±) = l(N t , F -N i>F )\t ± ) 

= ^AVlt±)=±~|t±)- (40) 

The occupation (or unoccupation) of this zero mode lead to N F = ±| charge and S z = \ 
bound to the defect^. Due to S^-conservation for the QAH state, the induced quantum 
spin number S z = | is also well defined. 

Thirdly, in the NI state, there is no zero mode on a single 7r-nux defect. As a result, all 
induced quantum numbers are zero. 

Therefore, in this part we find that the zero mode's on a 7r-fTux may also be used to dis- 
tinguish different quantum phases. These induced quantum numbers are shown in Table. 4. 





QSH 


QAH 


NI 


Zero mode's number 


2 


1 





N F 





i 

2 





S z 


i 

2 


1 

4 
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Finally we give a summary. When D ^ 0, because S z is a good quantum number, we 
can use the spin Chern number Q s to characterize different quantum states as shown in the 
Table.5. 





QSH 


QAH 


NI 


Spin Chern number Q s 


2 


1 





Edge state's number n s 


2 


1 





Zero mode's number 


2 


1 






For the TQPTs from one quantum state to another, the spin Chern number Q s will jump. 

C. Universal critical behavior of TQPTs 

Another feature of the TQPTs is the non-analyticity of the ground-state energy. The 
ground-state energy is defined as 

S f 

E = - ^2(e+(k) + £-(*)) = / + £ -(k))d 2 k, (41) 

k BZ 

where, BZ denotes the first Brillouin zone and 

'*» = ^ 2 ^ + ^*. + (f + *(«»*. + «»»M±^- (42) 
Here S is the area of the system. As an example, we show the TQPT of the model along 
the line with fixing D = 0.5 in Fig.l. As illustrated in Fig. 7, the ground state energy and 
its first and second derivatives are continuous for arbitrary t, while its third derivative is 
non-analytic at points t = 0.125 and t = 0.375, corresponding to phase transitions NI - QAH 
and QAH - QSH. It means that the TQPTs are third order. Similar third order TQPTs 
have been pointed out in other systems where the nodal fermion appears^. 

Let us explain why TQPTs are always third order. The energy dispersions near the 
quantum phase transitions are shown in Fig.8 and Fig. 9. On the line with fixing D = 0.5, 
one can see that there is one nodal point with zero energy in the Brillouin zone, (n, n), 
at the quantum phase transition between QSH state and QAH state. Thus the TQPT is 
dominated by nodal Dirac fermionic excitations at (tt, it). For the TQPT between QAH 
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FIG. 7: The ground-state energy and its first, second and third derivatives with respect to t, where 
D = 0.5 (which corresponds to the horizontal dash dot line in Fig.l). It is clear that E, ^ and 
are continuous functions, but is discontinuous at the transition points t = 0.125 and 0.375. 

state and NI state on the line with fixing D = 0.5, the nodal fermion is also pinned at the 
point (tt, 7r). 

Furthermore, near the critical points of the TQPTs, we get the low energy effective 
Hamiltonian near the point (tc, it), 

#i it m t r z + k y T x + k x r y (43) 

and 

#14 -> m X T z + k y r x - k x r y (44) 

where = —1 + At — D and = — 1 + At + D are the masses of the electrons with up-spin 
and down-spin, respectively. One can see that crossing the critical point D = 0.5, t = 0.125, 
the mass of the electrons with down-spin changes sign, -p^r — > — r~V Consequently, the 
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FIG. 8: The energy dispersion near the point of phase transition between QSH and QAH. (k x , 
k y : (-2vr, 2vr)) 




FIG. 9: The energy dispersion near the point of phase transition between NI and QAH. (k x , 
k y : (-2-7T, 2vr)) 

Chern number of the electrons with down-spin Q± changes from to 1 that corresponds to 
the TQPTs from NI state (Qt = Qi — 0) to QAH state (Q-j- — 0, Q± — 1). Similarly, crossing 
the critical point D = 0.5, t = 0.375, the mass of the electrons with up-spin changes sign, 
t— ! r, that leads to the Chern number of the electrons with up-spin changes 



T— T -> " 









from to 1, corresponding to the TQPTs from QAH state (Q^ = 0, Qi = 1) to QSH state 
(<5t = 1, Q4, = 1). For a special TQPT at (D = 0), the mass of the electrons with up-spin 
m-|- = — 1 + 4t is equal to the mass of electrons with down-spin = — 1 + At. The changes of 
mass sign will lead to the jumps of spin Chern number, AQ^ = AQ± = 1, that corresponds 
to a direct topological quantum phase transition from NI state (Q^ = Qi = 0) to QSH state 
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(<3f = 1, = 1)22. Fig. 10 is a scheme to illustrate the relationship between the jump of 
the spin Chern number and the change of mass sign. 



QAH (Q T = 0, Q l = 1) QSH (Q T = l, Q, = l) 




p bard (or up-spin electrons 

u band tor up-spin electrons 



m t > m t = m t < 

m, <0 m i <0 m, <0 



FIG. 10: Illustration of the relationship between the jump of the spin Chern number and the 
change of mass sign. 

On the other hand, due to the existence of the nodal fermion, the TQPTs with mass sign 
changes (77^4, — > — m-^) are always third order. One may use the low energy approximation, 
E{k) = \Jk 2 + rri^p and then get the ground state energy as 

E(m tti ) ~ -|m t J 3 + const t 4 (45) 

where const^ is a constant--. It is obvious that the third order derivative of E(jn^^) to 
m^ j, is discontinuous at the point m-j-^ = 0. 

In brief, for the TQPTs, NI - QAH and QAH - QSH, the changes of the mass sign 
— > —TWf^ reflect the jumps of the (spin) Chern numbers (AQ S 7^ 0). As a result, third 
order phase transitions is a universal feature of the TQPTs. 
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FIG. 11: The free energy and its first, second and third derivatives with respect to t, where D = 0.5 
(which corresponds to the horizontal dash dot line in Fig.l). It is clear that F, and 
are all analytic. 

D. Finite temperature properties 

At finite temperature (T / 0), the free energy F is denned by 

F = -T j -0^2 m(cosh(£ +(fc) +e_ (fc) )). (46) 

Here we set ks — 1. We calculate the derivatives of the free energy which are shown in 
Fig. 11. One can see that the free energy and its first, second, third derivatives are all 
analytic without singularity at T ^ 0. The results means that at finite temperature, there 
is no real phase transition, instead one gets crossovers. This result is consistent with that 
of the TQPT between QAH state and NI state in Ref.— . 

E. Stability of the TQPTs 

Finally we discuss the stability of the TQPTs. At the critical points, there always exist 
nodal fermion with zero energy at (tt, it). Thus at the critical points, we linearize the 
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dispersion at (tt, tt) and get a low energy effective Dirac Hamiltonian of nodal fermions, 

#1 = E *lt( r ^ + ^)**.t + E - ( 47 ) 
it /c 

it 

where the 7 matrices are defined by 7^ = r 2 and 7 y = r 1 . 

Then we consider effect of short range interactions. For example, one may add an on-site 
four-fermi interaction H int = U(^^f) 2 to Hi. Now we get an effective three dimensional 
Gross-Neveu model with the Lagrangian^ 

L GN = i¥ llx d^ + U(¥^) 2 . (48) 

It is known in large-N limit, the Callan-Symanzik function of (3(U) is 0(U) = + U. For 
a small interaction, U — > 0, the four-fermi interaction is irrelevant. That means the TQPTs 
is stable against the small on-site four-fermi interaction. Similarly, considering other types 
of short range interaction with S^- conservation, one may get the same results: they are 
irrelevant and the nodal fermions are stable. 

IV. TQPTS WITH T- SYMMETRY : D = 0, R + 
A. Global phase diagram 

In this section, to learn the TQPT of topological insulators with T-symmetry, D = 0, 
R 7^ 0, we focus on H 2 = H + Hr and show the effect of the Rashba term. 

From the global phase diagram in Fig. 12, one can see that there exist three quantum 
phases : QSH state, NI state and gapless semi-metal. The solid line divides the gapped 
phases and the gapless semi-metal. As shown in Fig. 13, when one fixes t at 0.375, the 
energy gap will decrease with the increase of R and close at R = 1. However, the phase 
boundary (the dashed line in Fig. 12) between the QSH state and NI state doesn't change 
by the Rashba term. 

B. Topological "order parameter" - Z2 topological invariant 

To see the topological properties of the Hamiltonian with T-symmetry but without S z - 
conservation, we calculate the edge states and the zero mode's number on a 7r-flux in different 
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FIG. 12: The global phase diagram at D = 0. 
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FIG. 13: The energy gap with respect to the Rashba term R at t = 0.375. 



quantum states. In the QSH state, from the numerical results, we find that there also exist 
two zero modes of a single tt-Aux defect. The charge density of a fluxon in the QSH state on 
a 24 x 24 lattice is shown in Fig. 14. In the NI state, there is no such zero mode on a single 
7r-ffux defect. On the other hand, the results of the edge states are shown in Fig. 15, from 
which we find that there are two edge states in QSH state, while in the NI state, there is no 
edge state. 

Above results indicate that the T-symmetry indeed protects the topological properties of 
the QSH state. Thus we calculate the Z 2 topological invariant (— 1) A = Ylt=i Ylm=i &m(ki) 
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FIG. 14: The charge density of a fluxon in the QSH phase on a 24 x 24 lattice with periodic 
boundary condition. The charge bound to the defect is ±e. 
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FIG. 15: Energy spectrum of the QSH state when the open boundary condition is imposed in 
x-direction. The parameters are R = 0.4, t = 0.5, fi = 2. 

in different gapped quantum state. For the four high-symmetry points, (0, 0), (0, it), (tt, 0), 
(7r, 7r), the eigenvalues below Fermi surface Ep of this system are 

l + 4t, 1, 1, l-4t, (49) 

respectively. TQPT occurs at t — 1/4. For t < 1/4, the four eigenvalues are all positive. 
Using the formula mentioned above, one can see that the Z 2 topological invariant (— 1) A is 
1 which means the system is in the NI phase. For t > 1/4, one eigenvalue is negative while 
the others are all positive. It is obvious that the Z 2 topological invariant (— 1) A is —1 which 
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means the system is in the QSH phase. 

Table. 6 shows the relationship between the Z 2 topological invariant (— 1) A and the topo- 
logical properties of the quantum states, 



D = 0, 


QSH 


NI 


Z2 topological invariant (— 1) A 


-1 


1 


Edge state's number 


2 





Zero mode's number on a 7r-fTux 


2 






C. Universal critical behavior of TQPT 



To study the universal critical behavior of the TQPT between QSH state and NI state 
without S^-conservation, we calculate the derivatives of the ground-state energy E. Here 
the ground-state energy becomes E = — J2k( £ +W + £ ~(k)) where 

e±(fc) = J sin 2 k y + (R^J (sin 2 k x + sin 2 k y ) ± sin 2 k y + + 2t(cos k x + cos k y )Y) 2 . (50) 

We show the TQPT along the line with fixing R = 0.2. As illustrated in Fig. 16, the third 
derivative is non-analytic at points t = 0.25, corresponding to quantum phase transitions 
NI - QSH. So the TQPT is also third order. 

The energy dispersion near the quantum phase transition is shown in Fig. 17. From the 
dispersion of the electrons near the critical points, we find that the TQPT is also dominated 
by the nodal fermion near the point (jr, it). 

Let us explain why the TQPT between QSH state and NI state are also third order. 
When R 7^ 0, there is no ^-conservation and one cannot define the spin Chern number. So 
we cannot use the jump of the Chern number to understand the TQPT. However, from the 
definition of the Z2 topological invariant, we find that near the TQPT, the mass signs of 
the two kinds of low energy fermions (m = — 1 + At) are determined by the Z2 topological 
invariant (— 1) A (See Eq H9|) . Consequently, for the TQPT between NI state ((— 1) A = 1) 
to QSH state ((— 1) A = —1), the sudden change of the Z 2 topological invariant will lead to 
the change of mass sign of low energy fermions. 

In brief, for the TQPT, NI - QSH, the change of mass sign come from the sudden change 
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FIG. 17: The energy dispersion near the quantum phase transition. The parameters are t = 0.25 
and R = 0.2. 
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of the Z 2 topological invariant. As a result, the TQPT between NI state and QSH state is 
always third order. 

V. TQPTS WITHOUT T-SYMMETRY AND ^-CONSERVATION : D ^ 0, R ^ 

In this section, we study the TQPTs of topological insulators without T-symmetry and 
^-conservation, D ^ 0, R ^ based on iJ 3 = H + H D + H R . 

2. . 1 1 1 . 1 . 1 
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FIG. 18: The global phase diagram at D = 0.5. 

We take D = 0.5 as an example (which corresponds to the horizontal dash dot line in 
Fig.l). The results are illustrated in Fig. 18. One can see that in the phase diagram there 
exist three phases : two NI states and one QAH state. The QAH state is robust while 
arbitrary small R will destroy the QSH state. To understand the disappearance of the QSH 
state, on the one hand the evolution of the energy levels of zero modes on a 7r-flux with R 
at the point D = 0.5, t = 0.75 is plotted in Fig.19. The black bars denote the continuum 
spectrum while the black triangles denote the energy levels of zero modes on a 7r-flux. One 
can see that the degeneracy of zero modes is lifted by the Rashba term. With increasing R, 
the energy gap of the continuum spectrum decreases, while the energy splitting of the zero 
modes AE increases. On the other hand, the evolution of edge states in QSH state with the 
Rashba term is shown in Fig. 20. When R ^ and D^0, due to the hybridization of two 
edge states, the edge states in the QSH state open an energy gap. In contrast, in the QAH 
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phase, the edge state and the zero modes are all stable. 




FIG. 19: The evolution of the energy spectrum of the system with a 7r-flux at D = 0.5, t = 0.75. 




FIG. 20: The evolution of edge states in QSH state with the Rashba term R at D = 0.5, t = 0.75. 
a, b, c and d correspond to R = 0, R = 0.2, R = 0.4 and R = 0.6 respectively. 

In this case (a system without T-symmetry and ^-conservation), the TKNN integer Q 
can be regarded as an "order parameter" to characterize different quantum states, Q = 
lBz(^ k x ^> z dk x dk y . Table. 7 shows the relationship between the TKNN integer Q and 
the topological properties of the quantum states, 
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QAH 


NI 


TKNN integer Q 


1 





Edge state's number 


1 





Zero mode's number on 7r-flux 


1 






Using similar approach in above sections, the universal critical behavior of the TQPTs are 
studied. The TQPTs are also dominated by nodal Dirac fermionic excitation at the point 
(ir, 7r). Since the TKNN integer Q (the Chern number) will jump crossing the TQPTs, 
AQ = 1, the mass sign of one low energy fermionic excitation changes. Consequently, the 
TQPTs are also third order. 

VI. CONCLUSION 

In this paper, based on a two-dimensional lattice model, we study the TQPTs between 
the QSH state, QAH state and normal band insulator and show their physical properties, 
including the edge state, the quantized (spin) Hall conductivity and the induced quantum 
number on a 7r-flux. There are common features of the TQPTs for different cases : the 
existence of nodal fermions at high symmetry points, the non-analytic third derivative of 
ground state energy and the jumps of the topological "order parameters". In particular, we 
find the symmetry protected nature of the TQPTs which is illustrated in Table. 8 : 



T-symmetry 


^-conservation 


" Order parameter" 


/ 


X 


Z 2 topological invariant (— 1) A 


X 


/ 


Spin Chern number Q s 


X 


X 


TKNN integer Q 



In Table. 8 the symbol / means the system is invariant of the given symmetry and x means 
not. From Table. 8, one can see that for a system with S^-conservation but without T- 
symmetry, the spin Chern number Q s plays the role of "order parameter" to characterize 
the topological insulators; for a system with T-symmetry but without S^-conservation, the 
Z 2 topological invariant (— 1) A becomes the "order parameter"; for the system without 
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S^-conservation and T-symmetry, it is the TKNN integer Q that becomes the "order pa- 
rameter" . 

In the end, we give a comment on the relationship between symmetry and topological 
invariants for the TQPTs. TQPTs are not be classified by symmetries, instead, they may 
be characterized by some topological invariants, such as the Chern number or Z2 topological 
invariant. However, the "symmetry" of the systems still plays important role : different 
topological quantum phase transitions are protected by different (global) symmetries and then 
described by different topological "order parameters". 

The authors thank Ying Ran and Xiao-Liang Qi for helpful discussions and comments. 
The authors acknowledges that this research is supported by NCET, NFSC Grant no. 
10874017. 
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